Liquid films with high surface modulus moving in tubes: dynamic wetting film and 

jumpy motion 
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We investigate the motion through a wet tube of transverse soap films, or lamellae, of high surface 
dilatationnal modulus. Combining local thickness and velocity measurements in the wetting film, 
we reveal a zone of several centimeters in length, the dynamic wetting film, which is significantly 
infiuenced by a moving lamella. The dependence of this influence length on lamella velocity and 
wetting film thickness provides a discrimination among several possible surfactant minimal models. 
A spectacular jumpy mode of unsteady motion of a lamella is also evidenced. 
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The motion of a liquid meniscus sliding over a wet 
wall is ubiquitous in many industrial and biological con- 
texts, among which dip coating [1], bubble or droplet 
motion in pores (e.g. in enhanced oil recovery) and mi- 
crofluidic channels [2], foam friction on solid boundaries 
[3], as well as lung diseases [4]. The problem has first 
been studied by Landau and Levich [5] and Derjaguin 
[g*] (LLD) for a soHd plate pulled out of a liquid bath 
at small velocity. For pure liquids, the dynamics is con- 
trolled by the dynamic meniscus that forms between the 
static meniscus and the wetting film withdrawn by the 
plate. The LLD model leads to a wetting film thickness 
^LLD = 0.945^cCa^^^, with the capillary length and 
Ca = r]U/^ the capillary number (7: surface tension, 
77: liquid viscosity, and U: plate velocity); the exten- 
sion of the dynamic meniscus along the plate scales like 
^LLD ^ ^cCa^/^. These predictions are based on the lu- 
brication approximation /i-LLD/^LLD ^ 1, and on a free 
shear boundary condition at the air/liquid interface. 

If the liquid phase is a solution of surfactants, the 
interfacial stress depends on the surface concentration 
of surfactants T and on the rheology of the surfactant 
layer. This stress reacts against the local area variation 
of the interface, with a dependence quantified by a sur- 
face dilatational modulus E. The friction against a wall 
of foams made of solutions with high E is well predicted 
by the limiting case of a locally incompressible interface 
[3] : the whole interface then moves with the foam, result- 
ing into a strong shearing of the wetting film between the 
foam and the wall. For single bubbles and foams, the to- 
tal interface area of each bubble is constant during the 
motion, which is compatible with the local incompress- 
ibility of the interface. 

In other geometries, like dip coating or motion of a 
lamella (a soap film across a tube), the total area of each 
connected part of the interface varies, and the location 
and extension of the zone where area variations occur 
becomes a central question. A classical assumption is 
that the area variations are localized in the static menis- 
cus, and that, in contrast with the model by Denkov et 



al. [3], the interface of the wetting film and of the dy- 
namic meniscus move with the plate. The hydrodynam- 
ical problem is then governed by the same equations as 
for pure Hquids, after a well adapted rescaling 0. Sev- 
eral models go beyond these limiting cases, by includ- 
ing surfactant adsorption/desorption [8], diffusion in the 
bulk Q or an intrinsic surface viscosity jlO)]. Recently, 
the full problem was numerically solved without the lu- 
brication approximation, which allowed to include the 
whole meniscus into the simulation, with surface diffu- 



sion 



Uj or surface diffusion and adsorption/desorption 



However, no clear experimental evidence allows to 
discriminate between these different models. 

Using a surfactant solution with high surface modulus 
and low bulk viscosity, we bring the first direct experi- 
mental evidence that a lamella moving in a wet tube can 
push the wetting film over centimetric distances, i. e. 
more than two orders of magnitude larger than I/lld- 
We thus introduce the concept of dynamic wetting film. 
We show that its lateral extension, that we call the influ- 
ence lengthy depends on the initial film thickness and on 
the lamella velocity. Moreover, we show that the various 
processes that may govern the influence length lead to 
very different scalings, discriminated by our experimen- 
tal data. Finally, we demonstrate a surprising jumpy 
behavior at high velocity, where the meniscus intermit- 
tently slips over its wetting film. 

Single lamellae are created and pushed at prescribed 
velocity /7 in a wetted vertical tube of inner diameter 
2a = 8.8 mm following the method of In order 

to maximize E^ we use the well characterized mixture 
of sodium lauryl dioxyethylene sulfate, cocoamidopropyl 
betaine, and myristic acid (MAc) in ultrapure water (so- 
lution Si), or the same with 40% wt glycerol (solution 
5*2), following the protocol of Q . The wetting film thick- 
ness h is measured by white light interferometry using the 
commercial spectrometer USB4000 (Ocean Optics). Ne- 
glecting multiple reflections within the film, the collected 
light intensity obeys /(A) = /o[q^ + /3 cos {Air nh / X)] with 
lo a reference intensity, n = 1.33 the optical index of the 



solution and A the wavelength. The film thickness /^fib(0 
is determined by computing the dominant Fourier com- 
ponent of 7(1/ A). Assuming a steady film profile in the 
lamella frame, we set h{x) = /ifib(^o — ^/U)^ with to the 
time at which the lamella is in front of the fiber. The 
direction x is oriented downstream, with x = at the 
lamella position (see FigUl right). A camera is synchro- 
nized with the spectrometer and records the film shape 
and velocity (Fig. [H left). 





FIG. 1: (left) Lamella passing in front of the fiber, (right) 
Scheme of the meniscus and notations used in the text. 



FIG. 2: Wetting film profiles obtained with different initial 
thicknesses and [/ = 0.6 cm/s (solution 5*1). The values ob- 
tained for Li and hi (i = 0, 1) after profile analysis are indi- 
cated resp. by thick segments and dashed lines. The central 
box corresponds to the meniscus and x = is the lamella 
position. Inset: Film profile obtained with U — 1.7 cm/s, 5*1 
(thick line), Us measured by tracer tracking for the same ex- 
periment, using a = 0.45 (dotted line) and Us deduced from 
the film profile using Eq. [2] (thin line). 



Typical wetting film profiles are given on Fig. [2l The 
signal modulation is destroyed for too high thickness gra- 
dients, hence the meniscus profile is not accessible. The 
wetting film thickness /lo, far ahead the lamella, depends 
on the waiting time since the previous lamella. Its evolu- 
tion by drainage is slow enough for to remains uniform 
and constant on the experimental space and time scales, 
in the range 1 to 20 /im. The thickness hi of the film 
deposited by the lamella does not need to be equal to 
h^. If ho hi, the meniscus volume changes, but slowly 
enough for the steady state assumption to remain valid. 

Most remarkably, the wetting film begins to swell sev- 
eral centimeters in front of the lamella, at a distance 
Lq. Some profile shapes are roughly exponential, but 
the thickest ones display a sharp transition between the 
flat wetting film and the swollen part of the film close 
to the meniscus, of typical thickness h ^ 2ho. The rear 
influence length Li is significantly smaller, but still centi- 
metric. We compute Lq and Li using a systematic fitting 
procedure leading to the characteristic distances marked 
on Fig. [21 Lq is plotted as a function of U and ho on 
Fig. [3l The data are somewhat scattered, but exhibit a 
clear tendency: Lq increases with ho, and decreases with 
both T] and U ; the best power law fit yields: 

(1) 



T ^ 7 l.OiO.lrr 

Lq (X ho U 



0.5±0.1^-1.2±0.1 



A similar law was obtained for Li (data not shown). 

In some experiments, the presence of tiny bubbles 
within the wetting film enabled the simultaneous mea- 
surement of the film thickness and of the surface velocity. 
These bubbles are convected at a velocity Ub = aug, with 
Us the surface velocity (the bar indicating velocities in the 
laboratory frame) and a < 1 an unknown friction param- 
eter. Comparison with the film thickness profile proves 
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FIG. 3: Influence length Lq as a function of the wetting film 
thickness ho for two bulk viscosities (full symbols: r]Si =1.0 
mPa-s, open symbols: 77^2 , = 4.0 mPa-s) and various velocities 
(in cm/s, o, •: 0.32, □, ■: 0.60, <>,♦: 1-15, >, ►: 2.35, <, ^: 
4.20). Inset: Lq as a function of Lo,d given by Eq. (|9|). 
Used parameter values: = = 0-1 N/m [3], D = 

10"^° mVs [3, Fo = 7 X 10"^ mol/m^ (from the area^ 
molecule ^ 25 [3l) and cq — cmac — 0.88 mol/m^ 
hence /ir ~ 8 /im. 



that, as the film begins to swell, the velocity rapidly in- 
creases from almost zero to a value of the order of U 
(Fig. [21 inset). The velocity gradient in the dynamic 
wetting film is thus close to U/h, and the friction force 
per unit length (in the z direction) can be estimated as 
fy ^ r]U{L()/ho -\- Li/hi). This estimate, in the range 
10 to 30 mN/m for all our experiments, is in good quali- 
tative agreement with the value extracted from the film 
images (Fig. ^ fy = 2joa/R ^^^ith 70 = 23.8 mN/m 
the equilibrium surface tension and R the radius of 
curvature of the lamella (data not shown). The friction 
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in the dynamic meniscus alone is given by the Bretherton 
law fy^B = 4.707oCa^/^ fv^ and is at most 1.6 mN/m in 
our range of velocities. It is thus negligible compared to 
the dynamic wetting film contribution for this solution. 

We also tested a SDS solution, with r] = 1.2 mPa-s 
and E < 1 mN/m and a solution of 12-hydroxy stearic 
acid with ethanolamine as a counterion [18], with £^ = 38 
mN/m and a rheothinning behavior with a high viscos- 
ity T] varying between 20 and 10 ~^ Pa.s for shear rates 
between 10 ~^ and 300 s~^. In both cases, no measurable 
influence length has been observed. This gives a strong 
hint that a large Lq is associated to solutions of high E 
and low rj. The influence length results from a compe- 
tition between the resistance of the surface against com- 
pression, and the viscous resistance of the bulk against 
shear. The former may arise from various microscopic 
mechanisms, governed by the surface viscosity or elastic- 
ity, coupled with the surfactant desorption rate or dif- 
fusion. We now investigate them in their simplest form, 
at the expense of a quantitative modeling; however, this 
scaling approach appears to be sufficient to identify only 
one scenario able to capture the main dependences exper- 
imentally observed, in the investigated parameter range. 

We model the region x > 0, ahead of the lamella, in the 
frame of the lamella. We consider a steady regime. Since 

<C a, the tube curvature is negligible and we assume 
an invariance in the z direction. In the dynamic wetting 
film but at a distance larger than Llld from the lamella, 
Marangoni effects dominate over capillary effects and we 
can neglect the Laplace pressure. The velocity u{x^ y) is 
thus oriented along x and varies linearly between — t/ at 
y = and Us{= Ug — U) at y = h. Mass conservation 
thus writes = dx[{us — U)h/2], hence: 



Us = U 1 



2ho 
h 



(2) 



in agreeement with experimental data (see Fig. [2linset). 

In the presence of the surface viscosity /i*, the tangen- 
tial stress continuity imposes, at = h{x) [19], 



(3) 



Ref. [14| shows that without MAc, the surface modulus is 
two orders of magnitude lower. Therefore, it is reasonable 
to assume that dxj is only related to dx^MAc^ hence to 
use a one-component model. Assuming for simplicity a 
linear relation for the surface tension 7 = 70 — £^(r — 
ro)/ro with To the equilibrium surface concentration in 
MAc and F = Fmac, we get the velocity: 



u{x,y) = - ( ii'dxxUs 

T] 



E 



dxr]y- U. 



(4) 



If the interfacial dynamics was dictated by surface vis- 
cosity rather than surface elasticity, Eq. (|4]) would be- 
come, ioT y = h and E = {]: 1-^Us/U = {ii*h/riU)dxxUs' 



From experiments, the variation of Ug is of the order of 
U (Fig. [21 inset), leading to the scaling: 



V 



1/2 



(5) 



which does not capture the observed dependence of Lq 
on the lamella velocity. 

We now neglect the interfacial viscosity to focus on the 
other contributions from the surfactants. The surfactant 
mass balance at the interface writes: 



dx{us^) = j, 



(6) 



with j the surfactant flux from the bulk to the inter- 
face. The dynamic wetting film is compressed by the 
moving lamella, inducing a typical surface concentra- 
tion variation AF > 0. Moreover, Eq. @ leads to: 

Us ^ -U- 



- M^dxT, from which we deduce u^dxT ~ U ^ 



and TdxUs ~ 
with ho ~ 10 



Ar Ehp 
r]Lo ' 



Experimentally, we get ^^^^ 



m, U 



1 cm/s and Lq ^ 3 cm, showing, 
as expected, that the convective term dx{usT) is domi- 
nated by the velocity variation, which tends to accumu- 
late surfactant in front of the lamella. The concentration 
gradient only reduces this effect. If we neglect this last 
term, we get the following scaling law for the convective 
term in (|6]), as already obtained in Q: 



dx{Tus 



EhoAT 



<0. 



(7) 



This term is balanced by the exchange term with the bulk 
j. It equals the diffusive flux in the subphase, and obeys 
a kinetic sorption law [19] that we assume linear. Then, 
j = —DdyC = k{c — T/hr) at = /i, with c the bulk con- 
centration in MAc, Co its equilibrium value, hr = Fq/cq, 
D the diffusion coefficient, and k the sorption velocity, 
both assumed to be the ones of the micelles in which 
MAc is solubilized [14]. If k is small, the exchange is 
limited by the desorption process, and c :^ cq. Then 
j ^ -kAT/hr, hence from (|6|) and (Q: 



Lo,a 



Ehohi 
Tjk 



1/2 



(8) 



This scaling is exactly the same as in the viscous case 
([5]), with an effective surface viscosity /i*^ = Ehr/k. 

If the desorption is fast, the surface is in equilibrium 
with the subphase and, close to the meniscus, the con- 
centration increase in the subphase scales as Ac{h) = 
AT /hr- Then, the flux directly depends on the bulk dif- 
fusive field: j = —DdyC. If the gradient is established 
over the diffusive distance hd ~ y^DLoJU, this leads to 
the scaling behavior: 



0,d ' 



Ehj 



2/3 



,2/3 



{DU) 



1/3- 



(9) 
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This is only valid if hd < Hq. Injecting (j9j) in 
the expression of hd^ we get the condition K = 
[EDhr/{r]U"^hQ)]^^^ < 1. With the experimental param- 
eter values (see caption of Fig. [3]), we get K ^ 2. If 
hd ^ ^0, the gradient dyC becomes small and j is esti- 
mated from surfactant conservation: j = dxij^ ucdy) o:^ 
dx{c Jq udy). From mass conservation, the flow rate 

f^udy is a constant, equal to —Uho. Hence, j ~ 
—UhoAc/ Lq^ and Lq in this convective regime is thus 
governed by: 



^0,c 



Ehr 



(10) 



which does not capture the observed dependence of Lq on 
ho. Overall, the only regime compatible with our experi- 
mental results is the diffusive one. Even if the exponents 
in ([T|) and (|9|) slightly differ, the rescaling by ([9]) shown 
in Fig. [3] (inset) provide a reasonable agreement with the 
experimental data, without adjustable parameter. More- 
over, we do not explain why Li is significantly shorter 
than Lq (Fig. [2]); this suggests that the interface resists 
more compression than extension. 

Finally, the most deformed lamellae deviate from 
steady motion by a striking scenario: they undergo pe- 
riodic "jumps" between phases of constant velocity (Fig. 
HI). These jumps are quick (less than 10 ms) and macro- 
scopic (of order 1 mm), which make them easily observ- 
able by naked eye. The curvature is partially released 
during jumps, and builds up again during the phases of 
constant velocity. Integrating Eq. (|3]) (with /i* = 0) 
along the front part of the dynamic wetting film yields 
7_ = 7o — /o^ rjdyudx ~ jo — ijULo/ho, with 7_ the value 
of the surface tension close to the meniscus. The largest 
values obtained for ULq/Hq are close to 70 , which shows 
that surface tension significantly decreases towards the 
meniscus. In this case, the linear relation between 7 and 
r does not hold anymore, which may explain the dis- 
crepancy between the experimental and theoretical ex- 
ponents. It is likely that the interface becomes unstable 
below a certain value of 7_ and collapses as in Lang- 
muir monolayers j2Q|, and that the jumps are macro- 
scopic manifestations of such an instability. This un- 
steady behavior is in marked contrast with that of SDS 
films, which never showed such jumps [13]. 

As a conclusion, our predictions of the length over 
which an air/liquid interface of a wetting film is entrained 
by a lamella may provide a quantitative criterion for the 
transition between tangentially immobile and mobile in- 
terfaces in foam/ wall friction With I a typical bubble 
size, these two limits correspond respectively to Lq ^ ^ 
and I/O <C ^, and the transition criterion Lq = £ may be 
expressed from our predictions, in terms of the material 
parameters of the surfactants. Experimental tests of this 
hypothesis are under way. 
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FIG. 4: Time evolution of the position of the meniscus (plain 
line) and of the leading point of the lamella (dashed line), in 
the laboratory frame and in the jumpy regime. Two snapshots 
show the shape of the lamella just before, and just after, a 
jump. 
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